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Abstract 

We describe the algorithm used to identify eharged traeks in the fixed-target charm- 
photoproduction experiment FOCUS. We begin by describing the new algorithm 
and contrast this approach with that used in our preceding experiment — E687. We 
next illustrate the algorithm's performance using physics signals. Finally we briefly 
describe some of the methods used to monitor the quantum efficiency and noise of 
the Cerenkov cells. 



3 



1 Introduction 



FOCUS is a fixed-target experiment concentrating on the photoproduction of 
charm that accumulated data at Fermilab from 1996-1997. It is a considerably 
upgraded version of a previous experiment, E687 [1]. In FOCUS, a forward 
multi-particle spectrometer is used to measure the interactions of high energy 
photons on a segmented BeO target. We obtained a sample of over 1 million 
fully reconstructed charm particles in the three major decay modes: 
K~n~^, i^~7r+7r^7r+, and K'-k^t:^ (and charge conjugates). We will 

refer to these as "golden modes" . 
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The FOCUS detector (see Figure 1) is a large aperture spectrometer with 
excellent vertexing and particle identification. A photon beam is derived from 
the bremsstrahlung of secondary electrons and positrons with an ^ 300 GeV 
endpoint energy produced from the 800 GeV/c Tevatron proton beam. The 
charged particles which emerge from the target are tracked by two systems 
of silicon microvertex detectors. The upstream system, consisting of 4 planes 
(two views in 2 stations), is interleaved with the experimental target, while 
the other system lies downstream of the target and consists of twelve planes of 
microstrips arranged in three views. These detectors provide high resolution 
separation of primary (production) and secondary (decay) vertices with an 
average proper time resolution of ~ 30 fs for 2-track vertices. The momentum 
of a charged particle is determined by measuring its defiections in two anal- 
ysis magnets of opposite polarity with five stations of multiwire proportional 
chambers. 
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Fig. 1. A schematic drawing of the FOCUS spectrometer. The target region consists 
of a segmented BeO target, the interleaved target silicon planes, trigger counters, 
and the 12 plane silicon tracking array. The spectrometer is approximately 32 meters 
long. 
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2 The FOCUS Cerenkov system 



Three multicell threshold Cerenkov counters are used to discriminate between 
electrons, pions, kaons, and protons. The Cerenkov system hardware was es- 
sentially unchanged from that used in E687 and is described in detail in Ref- 
erence [1] . A brief description of the Cerenkov system follows. 

There are three multicell thres hold detectors in the experiment, referred to as 
CI, C2, and C3. The detectors are operated at atmospheric pressure and in the 
threshold mode. The gases are chosen so that different indices of refraction (i.e. 
different light velocities) establish different momenta in which pions, kaons, 
and protons will begin to radiate Cerenkov hght (see Table 1). For our system 
the three pion thresholds were chosen to be 4.5, 8.4, and 17.4 GeV/c by use 
of appropriate gas mixtures. The photoclcctron yield ranged from roughly 2.5 
to 20 depending on the phototube and Cerenkov counter. 

The detector CI is the most upstream of the three Cerenkov counters, lying 
just beyond the first analysis magnet, between the first two PWC's (multi- 
wire proportional chambers) PO and PI. The gas used was a helium-nitrogen 
mixture, and the total length of the gas volume along the beam direction is 
180 centimeters. The Cerenkov detector C2 has the lowest threshold of the 
three detectors with a pion threshold of 4.5 GeV/c. The gas was pure N2O, 
and the total length of the counter gas volume along the beam direction is 
188 centimeters. The detector is located between PI and P2. 

The CI and C2 Cerenkov counters can detect all charged tracks that are gen- 
erally reconstructible in FOCUS. The 3rd Cerenkov detector C3 is located 
downstream of the second analysis magnet. Only higher momentum tracks 
make it through the aperature of this magnet, so C3 only helps in the identi- 
fication of these tracks. The counter is a helium threshold counter which was 
704 centimeters in length. 
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Table 1 



Cerenkov counter specifications. The momentum threshold for the three most rele- 
vant charged particles and the Cerenkov cone radius for a /? = 1 track at the image 
plane are given for each counter. 



counter 
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Cerenkov Radius 
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(cm) 
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15.9 


30.2 
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5.8 


CI 
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8.4 
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56.5 


90 


2.5 - 3.6 
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C3 


He 


17.4 


61.5 


117 


100 


9 


5.6 



3 The old versus new Cerenkov algorithm 



While the Cerenkov hardware used in FOCUS was essentially the same as 
E687, a completely new Cerenkov algorithm was written for FOCUS. This new 
algorithm will be referred to throughout this article by the acronym CITADL 
(for Cerenkov Identification of Tracks by an Algorithm using Digital Likeli- 
hood). Before describing the new algorithm, we briefly describe the previous 
algorithm known as LOGIC. For a more complete description of this algorithm 
see reference [1]. 

Unlike CITADL, whose decision is based on the individual firing pattern of 
all 300 cells comprising the FOCUS/E687 Cerenkov system, LOGIC based its 
identification on the overall firing status of CI, C2, and C3. LOGIC rendered 
a single identification indicating whether or not the track was consistent with 
the electron, pion, kaon, and proton hypothesis.f^ This decision was based on 
the track momentum and the Cerenkov light observed in the three threshold 
Cerenkov counters. A counter was declared "on" if any of the cells within the 
track's Cerenkov cone fired. A counter was declared "off" if no cells within 
the cone fired and a minimum number of expected photoelectrons (typically 
2.5) was expected under the pion hypothesis.P^ Otherwise the firing status 

Muons can only be effectively separated from pions over a narrow momentum 
range just below each counter's pion threshold. Both E687 and FOCUS had a sep- 
arate muon detection system to provide high quality muon identification. 

In order to save time, LOGIC computed the expected number of photoelectrons 
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for that counter was declared unknown and its information was removed from 
the final decision. The observed on or off firing status was then compared to 
whether or not the counter should have fired under a given hypothesis. This 
prediction was based solely on whether or not the track momentum exceeded 
an "effective" momentum threshold for that hypothesis.f^ 

Although the LOGIC algorithm was very effective at helping to isolate charm 
particles in E687, it did have shortcomings. LOGIC tended to discriminate 
against pions when one required positive kaon and proton identification. Much 
of logic's tendency towards light particle identification was intended given 
the goal of strongly suppressing pion backgrounds to the kaons found in 
Cabibbo favored charm final states. For example, any cell firing within the 
Cerenkov cone sufficed to declare a counter on. But if no cells fired, a sig- 
nificant amount of predicted light was required before that counter would be 
declared off. 

An unintended bias was due to accidental firings of the Cerenkov cells due 
to "noise." The noise was due to RF noise on cables, tube noise, and light 
from untracked, charged particles such as electromagnetic spray and photon 
conversions produced in the very intense photon beam. The electromagnetic 
noise source could be very serious for Cerenkov cells located in the center of the 
system where occupancies sometimes approached 25-50%. Both effects tended 
to assign Cerenkov light to tracks making them inconsistent with "heavy" 
particles such as kaons and protons. 

logic's tendency towards light particle identification both reduced the effi- 
ciency for kaon identification in Cabibbo favored decays and increased back- 
grounds for rarer Cabibbo suppressed decays such as D'^ vr+vr" or — > 
n^fi^u. In order to suppress the copious backgrounds from —>■ K^n^ or 

for each cell under the pion hypothesis unless the track was under the pion threshold 
for the Cerenkov counter. If a track was below pion threshold, the light yield was 
computed under the electron hypothesis. This allowed the Cerenkov system to help 
in the identification of electrons. 

The "effective" threshold was slightly higher (~ 10%) than the actual threshold 
in order to crudely take into account the gradual rise in the expected light yield 
with momentum above threshold. 
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K fi^u, one would typically require that the pion had a Cerenkov 
response which was inconsistent with that for a kaon. 

While studying Cabibbo suppressed states, we used D^^s skimmed from a 
sample of D*^ D^tt^ with no Cerenkov cuts to measure the fraction of kaons 
which passed the pion cuts. Typically 5% of kaons were misidentified as pions 
in E687 by the LOGIC Cerenkov algorithm. Because of the inflexibility of the 
LOGIC algorithm, one would need to redesign the internal cuts to minimize 
the misidentification of kaons as pions and re-run the algorithm from tapes 
that had the required Cerenkov ADC information. Although, in principle, 
LOGIC could be re-run with other internal cuts, that was not a practical 
option since data summary tapes typically did not contain the Cerenkov ADC 
information. 

CITADL is primarily motivated by the desire to produce a more flexible 
Cerenkov identification algorithm than LOGIC. In fact, the overall perfor- 
mance of CITADL was significantly better than that of LOGIC, primarily 
because CITADL allows for the possibility of accidental firing of Cerenkov 
cells. Rather than making a hard decision, on whether or not a track was 
consistent with a given hypothesis, CITADL returned relative likelihoods that 
the track had a Cerenkov pattern similar to that expected for the electron, 
pion, kaon, or proton hypothesis. One could then, for example, put a mini- 
mum cut on the likelihood ratio that the kaon hypothesis is favored over the 
pion hypothesis in order get sufficiently clean kaons to do the required physics. 
Unlike LOGIC, very few cuts were required to be "hardwired" in the CITADL 
algorithm. 

Like LOGIC, CITADL only uses the on/off status of Cerenkov cells rather than 
their pulse height in identifying particles. This decision made the computation 
of likelihoods simple since a cell's firing probability is given by the Poisson 
probability (1 — exp(— /x)) where is the expected number of photoelectrons 
under the given particle hypothesis.p^ 

This assumes that the gains and thresholds are set such that a single photoelectron 
wiU produce an ADC count in excess of the threshold required to caU a cell on. 
Under this assumption a cell will fire unless photoelectrons are observed when fj, 
are expected. The Poisson probability of getting zero photoelectrons is exp(— ^u). 
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CITADL constructed a log likelihood variable based on the firing probability 
for all Cerenkov cells that a given track could potentially affect all cells within 
the track's P — 1 Cerenkov cone. Assuming for the moment, that a cell only 
fired in response to Cerenkov light, if the cell fired, and // photoelectrons 
were expected, the log likelihood was incremented by log(l — exp(— //)); while 
if the cell failed to fire the log hkelihood was incremented by log(exp(— //)). 
Cells which were inside more than one track's Cerenkov cone were considered 
"confused" and excluded from the sum. The hkelihood returned by CITADL 
is similar in spirit to the traditional continuous likelihood used in fitting. The 
only difference is that each event has only two outcomes — on or off. For this 
reason, we call it a "digital" hkelihood. 

CITADL returns its identification in the form of like variables which we will 
call We , W^, Wk, and Wp. They are defined by Wi = -2 E^'"' log Pj where Pj 
is the probability for the observed outcome (on or off) for the j'th cell under 
each of the 4 particle hypotheses. One would typically require that potential 
charm decay kaons pass a minimum cut on a likelihood difference variable 
such as AWk = W^^ — Wk- A large AWk implies that the kaon hypothesis 
is significantly favored over the pion hypothesis. Unlike the case in LOGIC, 
there is no need to introduce "effective" thresholds , since the /i dependence on 
momentum is explicitly taken into account. There is also no need to declare a 
minimum number of photoelectrons required for a Cerenkov decision. If a very 
small number of photoelectrons discriminated the two hypotheses, CITADL 
returns likehhood differences close to zero. 

In computing the log likelihood, we take into account the probability that a 
given Cerenkov cell fires accidently due to noise as well as firing due to a given 
track. We determined the accidental firing rate by measuring the fraction of 
times a Cerenkov cell would fire, even if it were outside of the P — I Cerenkov 
cone of all observed tracks. A typical plot of the accidental rate as a function 
of cell number for one of the runs is shown in Figure 2. The accidental rate 
varied considerably and for central cells was very large. It is very easy to 
incorporate accidental firing rates in the firing probability. The prescription is 
-Pfire = a + (1 — exp(— yu)) — a (1 — exp(— //)) where a and // are the accidental 
rate and the number of photoelectrons expected for the given cell. We found 
that a was often proportional to the beam intensity — especially for cells near 
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the beam axis. CITADL included this effect as well. The inclusion of realistic 
accidental rates significantly improved the performance of our new algorithm 
relative to LOGIC. 
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Fig. 2. The fraction of times (in percent) that a cell in C2 fires when no detected 
track's /3 = 1 Cerenkov cone impinges on the ceh. These data were accumulated over 
a single run. Although most of the cells have an accidental rate of a few percent, 
cells located near the beam axis have accidental rates as high as 40%. 



4 CITADL Performance 



The very high statistics FOCUS data set provided numerous checks of the 
performance of the Cerenkov system and the CITADL algorithm - often on a 
run-by-run basis. These checks used decays into final states of known daugh- 
ters. The decays Ks tt+tt^ provided a very pure and highly copious source 
of pions, consisting of 15,000 decays in each of our nearly 6000 runs. This 
sample was large enough to provide an accurate photoelectron re-calibration 
for nearly all of the 300 cells in the Cerenkov system. 

Although not nearly as copious as our Ks sample, the decay A pn^ provided 
a clean sample of proton and low momentum pion decays. Finally the decay 

Reference [2] describes the method used to reconstruct the Ks and A topologies 



11 



K^K^ was used to measure the Cerenkov identification of kaons on a 
run-by-run basis.[^ Tlie run-by-run fraction of misidentified daughters from 
the -ft's'jA, and decays for several Cerenkov cuts was used as a stabihty 
monitor of the Cerenkov system. Figures 3 and 4 show examples of these 
"misidentification" monitors. 

Misidentification of pions from Kg ^ tt tt 



'^%0^''^^i^Km'i>'mmti6^'m^^<i> MiniiiMtiiiifiiftfiintr ! 

" .... I .... I .... I .... I . 

50 8000 10000 12000 14000 

Run number 

Fig. 3. The fraction of times that a pion from Kg ■k^tt~ is misidentified as a 
kaon, proton, or electron for three different CITADL cuts. Each point is averaged 
over 25 runs. 

We also found that it was possible to use golden mode charm as a monitor 
of Cerenkov performance. Figure 5 shows a 405,000 event golden mode charm 
sample obtained (using about 75% of our data) without any Cerenkov cuts. 
A selection of cuts on vertex detachment, isolation, the D*^ — mass differ- 
ence, and momentum were used to obtain this reasonably clean sample. Also 
shown are sideband regions used for background subtraction. Figure 6 shows 
the likelihood difference /^Wk = W^, — Wk for the kaon and pion daughters 
from these background subtracted charm decays for tracks with two ranges 
of momentum. For convenience, we will call the variable AWk = — Wk 
"kaonicity". A positive kaonicity implies that a given track is more likely to 
be a kaon as opposed to a pion. 

in FOCUS. 

To obtain a clean enough cp sample to make a meaningful background subtraction, 
we required that one of the two kaons was Cerenkov identified. 
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Misidentification of protons from A ^ P tt 
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Fig. 4. The fraction of times that a proton from A ^ pvr is misidentified as a light 
particle for three different CITADL cuts. Each point is averaged over 25 runs. 
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Fig. 5. Invariant mass plot for the three golden mode decays 

K^TT^ , K^n^TT^TT^ , and — s- K~'K~^7r~^. The reconstructed 
mass was shifted by 5 MeV/c^ so that its peak will reconstruct in the same place 
as the peak of the D^. This data has vertex quality and kinematic cuts only. No 
Cerenkov cuts were used. The vertical lines denote signal and sideband regions 
which will be used to make a background subtraction. 

Figure 6(a) shows the kaonicity distribution for charm kaons and pions in a mo- 
mentum range above the pion threshold of C2 (the lowest threshold counter) 
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Fig. 6. The log likelihood difference Wt^ — Wk distribution obtained from back- 
ground subtracted kaons (x's) and pions (diamonds) from the golden mode charm 
signal shown in Figure 5. The pion distributions were rescaled to have the same 
area as the kaon distributions. Fig. (a) Tracks with momenta in the range 5 < P < 
60 GeV/c. Fig. (b) Tracks with momenta in the range 9 < P < 16 GeV/c. There 
are off scale spikes in the 0-bin consisting of 20,000 and 4,500 events for Fig. (a) 
and (b). 

but below the kaon threshold of C3 (the highest threshold counter). Outside 
of this momentum range, the FOCUS Cerenkov system is incapable of much 
K-TT separation and the kaonicity distribution is strongly peaked near zero.p^ 
Figure 6(b) shows the kaonicity distribution in the more restricted range from 
9 to 16 GeV/c. In this range kaon-pion discrimination is particularly effective 
since it lies above the pion threshold for CI but below the kaon threshold of 
C2. 



Figure 6 shows that, even though the likelihoods are constructed from the dis- 

CITADL offers some slight K-vr discrimination outside of this range since it can ex- 
ploit the momentum dependence of photoelectron yield beyond the C3 kaon thresh- 
old: i.e. the threshold is not infinitely sharp. 



14 



Crete firings of Cerenkov cells, the kaonicity distribution for kaons is reasonably 
continuous except near AWk = 0. As Figure 6(a) shows, averaged over the 
accepted charm momentum spectrum, pion backgrounds to kaons can be very 
effectively eliminated while still maintaining high efficiency for charm kaons. A 
cut just above kaonicity of zero rejects a large fraction of pions. The fraction of 
background pions dies away exponentially with the kaonicity cut above zero. 
Over the more restricted range from 9 to 16 GeV/c, where cells from both CI 
and C2 discriminate pions from kaons, the AWk distribution shows a sigifi- 
cantly larger average kaonicity. One can make a very stringent kaonicity cut 
to suppress pion backgrounds and still maintain good efficiency for real kaons. 

The situation for pion identification is essentially the mirror image of that 
for kaons. The contamination of kaons into the AWk < region falls off 
exponentially in AWk, while the pion spectrum extends below AWk < —20. 
In the region from 9 to 16 GeV/c, where both CI and C2 discriminate pions 
from kaons, the average kaonicity of pions becomes significantly more negative 
permitting one to make more stringent cuts to reduce misidentification. 

4-1 Understanding the kaonicity distributions 

We begin by discussing the kaonicity distribution for kaons in Figure 6. As 
described below, we use a simplifiedp^ model to conclude that the positive 
half of the kaonicity distribution is controlled by twice the total number of 
photoelectrons in those Cerenkov cells which discriminate pions from kaons 
and the negative half of the spectrum depends on the accidental firing rate 
but is damped exponentially in kaonicity. 

When a track is assigned a non-zero kaonicity, there are some Cerenkov cells 
that would be expected to fire if the track were a pion and would be expected 

One simplification is that the model assumes that all pion-kaon discrimination in 
CITADL is due to cells which have pion thresholds below the track's momentum but 
kaon thresholds above the track's momentum. In fact there is extra discrimination 
for cells with kaon thresholds below the track momentum since the probability that 
the cell will fire under the kaon hypothesis is less than the firing probability under 
the pion hypothesis. 
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to fail to fire if the track were a kaon. If AWk > 0, these "discriminating" 
cells did not fire. Let /i be the sum of the photoelectrons expected for pions for 
these "discriminating" cells. Assuming a small accidental rate, the probability 
that these cells would fail to fire under the kaon hypothesis is essentially 1, 
and under the pion hypothesis is exp(— /i). The kaonicity is defined in terms 
of the probability that the light pattern agrees with the kaon hypothesis di- 
vided by the probability that the pattern agrees with the pion hypothesis or 
AWk = —2 log (PttZ-Px)- After inserting these non- firing probabilities, we have 
AWk = 2/i. According to this model, the maximum kaonicity of ~ 30 means 
that the FOCUS Cerenkov system provides at most 15 photoelectrons which 
discriminate between kaons and pions. In the momentum region between 9 and 
16 GeV/c, Figure 6(b) shows that the yield of discriminating photoelectrons 
is much larger than over the full range from 5 GeV/c to 60 GeV/c. This is 
because pions with momenta in this range should fire the cells of both CI and 
C2 whereas kaons should not. The photoelectron yield for (3 = 1 tracks in C2 
is typically in excess of 11 — larger than that for CI or C3. 

If CITADL assigns a track AWk < 0, there must be discriminating cells which 
fired making the pion hypothesis more likely than the kaon hypothesis. For 
real kaons, such as those displayed in Figure 6, this can only happen due to 
accidental firing. Denote the probability that noise fires a discriminating cell 
by a. In the limit where there is a reasonable number of discriminating photo- 
electrons, the probability that the cells will fire under the pion hypothesis will 
approach 1. Hence if a cell accidently fires for a kaon, CITADL will report a 
kaonicity of AWk = —2 log^P-^/PK) = 2 log(a) where a is the accidental firing 
rate. The probability that a kaon track will actually fire a discriminating cell 
will of course be a = exp^AWKf^)- One therefore expects the roughly expo- 
nential fall off in negative kaonicity for real kaons which is observed in Figure 
6.p^ To summarize, a positive kaonicity value is essentially twice the number 
of photoelectrons which discriminate kaons from pions at the momentum of 
the kaon; while the negative half depends on the distribution of accidental 



The distribution of negative kaonicities for kaons also depends on the distribution 
of noise rates for cells in the Cerenkov system. Deviations from an exp(AVl^i^/2) dis- 
tribution are therefore expected since the distribution of accidentals is nonuniform 
as shown in Figure 2. 
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firing rates and is suppressed by a factor of exp{AWK/2). 

We next turn to a discussion of the AWk distribution obtained for the pi- 
ous shown in Figure 6. For negative kaonicities {AWk < 0), the kaonicity 
distribution exhibits considerable structure but when Al^;<: > it dies ex- 
ponentially with kaonicity. The exponential fall-off in the Al^/^ > region 
is due to discriminating cells not firing for the pion thus causing CITADL to 
prefer the kaon hypothesis. For pions, this should happen with a probability 
of — exp(— /x) while for kaons it will occur with — 1. Such tracks will 
therefore be assigned a kaonicity of AWk — — 2 log(P7r/-Pft') — We thus 
expect a kaonicity distribution for AWk > given by the product of the spec- 
trum of discriminating photoelectrons times the probability of the pion not 
firing the discriminating cells. This leads to a nearly exponential distribution 
since the probability of a pion not firing the discriminating cells is given by: 
exp(— ^) = exp(— In fact, the kaonicity spectrum of pions in the 
region 0.5 < AWk < 10 is well fit to the form exp(-0.4 AWk)- 

We next consider the AWk < half of the kaonicity spectrum for pions. 
In this region, the pion hypothesis is favored over the kaon hypothesis if the 
discriminating cells fired. Assuming a relatively large number of discriminating 
photoelectrons, the probability the cells will fire is close to P-^- = 1. Under 
the kaon hypothesis the cells would only fire due to accidentals which would 
occur with a probability of Pk = a. The pion would then be assigned a 
kaonicity of AWk = — 2 log(P7r/Pft:) = 21og(a). If several discriminating cells 
fire, the kaonicity distribution will be incremented by several multiples of 
2 log(a). Indeed several peaks are present in the AWk < spectrum of Figure 
6(b) which appear at multiples of approximately 6.25 which implies a typical 
accidental rate of a = exp(— 6.25/2) — 0.044. This estimated a is consistent 
with the typical accidental rate shown in Figure 2. 

To summarize, the AVF^- > spectrum for pions is controlled by the number 
of discriminating photoelectrons and is damped by a Poisson inefficiency factor 
of exp(— Al^/2) while the AWk < region is controlled by multiples of twice 
the log of the accidental firing rate. 

This model, among other things, explains why the pion-kaon separation in the 
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momentum range of Figure 6 (b) is so much better than over the complete 
momentum spectrum. Because CI and C2 both discriminate in this narrow 
momentum range, there are more Cerenkov cells available to discriminate be- 
tween the pion-kaon hypothesis. Hence there are more discriminating photo- 
electrons which increases the average kaonicity for kaons, and more multiples 
of —2 log(a) which further decreases the (negative) kaonicity distribution for 
pions. 

4-2 Using Cerenkov Information to Reduce Charm Backgrounds 

As Figure 5 shows, it was indeed possible to get reasonably clean charm signals 
without the use of Cerenkov information. However, many FOCUS analyses em- 
ployed Cerenkov cuts as an effective way of increasing signal to noise, while 
maintaining reasonable efficiency. Figure 7 illustrates the effectiveness of kaon 
and pion Cerenkov cuts for K^-k^-k^ti^ events selected by requiring 

that the secondary to primary vertex detachment exceeded 9 standard devi- 
ations. No Cerenkov cuts were used in the initial selection. The kaon cut is 
on "kaonicity" or the log likelihood difference AWk = W^, — Wk discussed 
previously. The pion cut is based on a pion consistency variable which we 
will call "piconicity" or APV^ = Wmin — W-,, , where PVmin is smallest negative 
log likelihood of the 4 particle hypotheses. The APV^ cut is placed on all D 
decay pions and is meant to ensure that no pion being considered as a charm 
daughter is grossly inconsistent with the pion hypothesis.p^ A cut such as 
AWjr > —2 means that none of the other 3 particle hypotheses is favored 
over the pion hypothesis by more than a factor of exp(2/2) = 2.71. For the 
D'^ — >■ K~7T~^7r~^7T~^ sample the requirement that AWk > preserves 84% of 
the yield while increasing the signal to noise by a factor of 6.2. The more 
stringent AWk > 2 and AW-,, > —2 preserves 75% of the uncut signal yield 
but increases the signal to noise by a factor of 16. 

27 We generally use a consistency cut rather than demanding that the pion is favored 
over both the kaon and electron hypothesis since the momentum range at which 
pions can can be distinguished from electrons is below 17 GeV/c for tracks traversing 
all three Cerenkov counters and below 8.5 GeV/c for 3 chamber pions which traverse 
only CI and C2. 
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Fig. 7. Illustration of the effectiveness of Cerenkov cuts in reducing backgrounds to 
Z)0 ^ K- TT+TT+vr . Note the logarithmic scale. The upper curve has no Cerenkov 
cuts. The second histogram requires AWk > 0. The third requires AWk > 2. The 
fourth histogram requires AWk > 2 and AWt^ > —2. A considerable improvement 
in the signal to noise is evident with only moderate loss in efficiency. The fitted 
signal yields in these plots are 15307, 12783, 11699, and 9996 respectively. 

One of the goals of the CITADL algorithm was to be much more efficient than 
LOGIC in suppressing the number of kaons which are misidentified as pions to 
enable us to more effectively study Cabibbo suppressed decays. An example of 
such a process is the tt^tt^ which is plagued by a large background from 

misidentified K^tc^ decays which occur with a branching ratio that is 

approximately 25 times larger than that of tt+tt^. Figure 8 compares 

the dipion mass spectrum from the published E687 signal to a version from 
half of the FOCUS data set. The E687 sample used LOGIC. The FOCUS 
sample required Wk — Wt, > 3 for both pions in order to significantly reduce 
K^-K^ contamination. Much of the improvement in event yield is due 
to the fact that FOCUS took roughly a factor of 15 times the E687 data set. 
We have also required that the dipion vertex be outside of the FOCUS target 
material to further increase our signal to noise relative to E687. However, the 
CITADL algorithm is responsible for the significant reduction in FOCUS data 
of the reflection from misidentified K~^tt~ relative to the tt^tt^ 

signal compared to what was achievable in E687. 
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Fig. 8. (Left) E687 signal for vr^vr from Reference [3]. The massive distor- 

tion in the background at lower masses is due to contamination from misidentified 
£)0 ^ K+TT-. (Right) FOCUS signal for vr+vr- from half of our data set. 

Figure 9 is an example of a plot used to gauge the effectiveness of a set of 
Cerenkov cuts on the pions and kaons from a very small sample of 
K~T{^Tx^ decays. The data satisfied our standard skim cuts for this mode: a 
good quality vertex intersection (CL > 1 %) , a kaonicity cut of /^Wk > 0.5, 
and a secondary to primary detachment exceeding 2.5 standard deviations 
{i/a > 2.5). This particular plot used the sample of decays which verti- 
cized outside of the target material and target microstrip system to remove 
backgrounds from multiple interactions. We show the yield versus signal to 
noise for 2 detachment cuts, and a sequence of Cerenkov cuts on the kaons 
and pions.P*] Figure 9 shows that the "piconicity cut" is essentially as effective 



Both the yields and signal to noise were based on fits to a Gaussian signal over a 
polynomial background. We define the signal to noise ratio as the ratio of the fitted 
number of signal events at the peak over the fitted number of background events at 
the peak mass. 
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a cut as the kaonicity cut. Figure 9 also shows that the Cerenkov cuts increase 
the signal to noise by a nearly constant factor at the two detachment cuts 
being considered. 
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Fig. 9. Illustration of the effectiveness of Cerenkov cuts in reducing backgrounds to 
— > K^TT^TT^ which verticize outside of the target microscrips and target mate- 
rial. We form a "cut tree" by plotting the signal yield versus S/N for two different 
detachment cuts and several cuts on kaonicity and pion consistency (piconicity). 
The kaonicity cuts (the main trunks) range from AWk > 0.5,1,2,3,4,5. The pi- 
conicity cuts (the branches) are "no cut", AWt^ > —10, —9, —3. Only 2 % of the 
our complete data was used for this plot. 



5 Calibration and Monitoring 

We made a large number of plots while reconstructing our data to monitor 
Cerenkov system performance. Examples of such plots which we have already 
discussed include Figures 2, 3, and 4 which serve as monitors of accidental 
rates, and Vee [Ks or A) daughter misidentification. We found that our most 
powerful calibration tool used pions from Ks tt+tt^ since roughly 15,000 
clean Ks^s were reconstructed during each of our ^ 6000 data runs. These 
pions were used to study how well we could predict the firing rate for cells. An 
example of such a study is shown in Figure 10 which plots the observed average 
firing rate versus the predicted. These particular plots are summed over all 
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cells in each Cerenkov counter. They represent the statistics on a single run. 
The prediction depends on the assumed cell geometry, the accidental firing 
rates, the relative quantum efficiency of each tube, and the validity of the 
analytic model used to quickly compute the fraction of Cerenkov light falling 
within the cell boundary. This model was known to have problems for cells 
on the planar mirror boundaries of CI and the planar mirror apex of C2. It is 
clear from Figure 10 that the light predictions were imperfect. We believe that 
the impact of these imperfections was slight on overall Cerenkov identification. 




0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 



P (predicted) p (predicted) 




P (predicted) 

Fig. 10. Plots of the actual firing rate as a function of the predicted firing rate for 
all cells within the given counter. We use pions from the decay Kg vr^vr". This 
set of plots was obtained in run 13955. 

We worked diligently to insure that CITADL used good information on the 
performance of each of the 300 cells comprising the FOCUS Cerenkov informa- 
tion. The two critical ingredients are the photoelectron yield for a, (3 = 1 track 
and the "accidental" firing rate. We found that the performance of some pho- 
tomultipliers varied significantly during the 12 month run. Occasionally this 
was due to the maintenance of the Cerenkov system, such as changing bases 
and phototubes. More often, small shifts in the ADC pedestals would create 
a dramatic apparent increase in the "noise" level of the tube, which could 
be easily corrected by raising the minimum ADC count required by CITADL 
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to call a cell "on." Although special calibration runs were taken in order to 
understand the Cerenkov system, we found that the best monitoring of the 
Cerenkov system was obtained through regular data taking. In this section, 
we describe some of these in situ calibration and monitoring methods. 



5.1 Photoelectron Calibration and Monitoring 



We developed a powerful way of monitoring the photoelectron yield as a func- 
tion of time for nearly all of the 300 cells comprising the FOCUS Cerenkov 
system. This cahbration method fit for the P — I photoelectron yield of each 
cell by minimizing CITADL likelihood Wt^ for pions from Ks — > n'^n'. The 
Wjr was incremented for each Ks pion daughter that was predicted to leave at 
least 0.1 photoelectron in a given cell. The hkelihood was incremented using 
background subtraction weight of 1 if the Kg mass was in the signal region 
and —1 if the mass fell in symmetrically placed, half width sidebands. Sep- 
arate Z^Wtt sums were computed for different assumed P — 1 photoelectron 
yields for the given cell which ranged from 20% to 160% of the nominal phot- 
electron yield. Typical log likelihood versus photelectron ratio plots obtained 
in a single run are shown in Figure 11. With the statistics available in a single 
run one could get an adequate re-calibration of the inner Cerenkov cells which 
are struck most often. 

As a stability monitor, we also summed the Wt^ from all of the Ks pions in 
a given run over all of the cells in our three Cerenkov counters. The minima 
of these grand likelihoods are plotted as a function of run number in Figure 
12 and give an overall scale factor for each of the three counters . Apart from 
some photoelectron fluctuations unique to CI in the early running, Figure 12 
shows that the average photoelectron yield from all three Cerenkov counters 
tended to fluctuate together in a way which we learned was correlated with 
changes in barometric pressure. 
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PE calibration for run 13955 
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Fig. 11. Calibration curves for three cells in the Cerenkov system. We plot J2^it 
for pion daughters from background subtracted Ks tt+tt" decays which strike 
the cell summed over all such pions in run 13955. This likelihood is computed for 
15 assumed photoelectron ratios relative to the /? = 1 yield used by CITADL. The 
example cells are in the inner section of each Cerenkov counter 

5.2 Noise Calibration and Monitoring 

CITADL makes direct use of the observed accidental firing rate for Cerenkov 
cells which is generally ^2 — 5% for most cells but can be quite high 40%) 
for cells close to the beam axis. We believe much of this noise for central cells is 
due to the high rate of e^e~ pairs that accompany our hadronic photoproduced 
events. Most of the very "noisy" cells showed an accidental rate which was 
roughly proportional to the instantaneous beam intensity. Figure 13 illustrates 
this point for two Cerenkov cells. A few words are in order. As a measure of 
the accidental rate, we use a variable called the "Poisson Accidental Rate". 
This "Poisson" rate (/la) is related to the accidental rate a (or fraction of times 
a cell fires when no track is pointing at it) via 1 — exp(— yU^) = a. We assumed 
in our treatment that fia rather than a was linear in the beam intensity.p^ 
We found that the average number of embedded pair tracks per triggered 

For small a , a but clearly as a ^ 1, a cannot continue to grow linearly with 
intensity. 



1600 
t= 1400 
] 1200 

1000 
BOO 



' r ' ' ' r ' ' ' r ' ' ' r ' ' ' I ' 
Cl-65 



. I .... I .... I .... I . 



BOO 
700 
600 
500 
400 



24 



1.4 
1.3 
1.2 



03 



1 



LU 0.9 

0.7 
0.6 









1 CI yM^-^t- •• 


V V ^"4 












• 








I . MfV^i^**^ ■'^ C3 






ISA-" . 



6000 



8000 10000 12000 

Run number 



14000 



Fig. 12. Plot of the relative photoelectron yield for all three counters as a function 
of run number. Each point is a 10 run average and the data for each counter are 
offset vertically for clarity. The photelectron yield relative to that assumed in the 
calibration is always within 20% of unity until run 9000 and within 10% thereafter. 
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Fig. 13. Plots of Poisson accidental rate versus instantaneous beam intensity for 
two Cerenkov cells. 

event formed a very convenient monitor of our average instantaneous beam 
intensity. Embedded pair tracks are from a Bethe-Heitler beam conversion 
which happen to lie within the resolving time of the chambers and microstrips 
when triggered by one of our event triggers. Pair tracks are easy to identify 
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since they are consistent, within their expected multiple scattering, with being 
produced along the beam axis.P^ 

We found that the average number of pair tracks in coincidence with an event 
trigger was a good measure of the instantaneous beam intensity which included 
the sometimes dramatic effects of spill non- uniformity. We used linear fits of 
the Poisson accidental rate versus (pair tracks), such as those in Figure 13, 
to model the noise response for each of the 300 Cerenkov cells. Often several 
parameterizations were made for a given cell to cover run dependent changes 
in the accidental rate. 

In fact the beam intensity varied significantly between the spills and even 
within a spill.[^ To optimize Cerenkov algorithm performance, we devised a 
method to estimate the beam intensity directly before each recorded event 
and use this intensity measurement to estimate the accidental rate for each 
cell. This method used several scalers which were recorded on our raw event 
tapes during each second level, event trigger. One of these scalers counted the 
accelerator RF clock, the others counted the hits in two scintillation counters 
used in the first level trigger. The scalers were effectively reset each time we 
read an event out. The time rate of change of either of these two scintillation 
counters formed a direct measurement of the beam intensity right before the 
actual event. Figure 14 shows that both the average number of embedded pair 
tracks, and the accidental Cerenkov cell firing rates are strongly correlated 
with the scaler derived beam intensity even within a single spill of Run 7185. 



6 Summary and Conclusions 

In this article we describe a likelihood based Cerenkov algorithm used to 
identify charged particles in the FOCUS charm photoproduction experiment. 

The average number of embedded pair tracks was « 1 per event throughout much 
of the FOCUS running owing to our very high intensity. 

Typicahy we had one spih of protons every minute which lasted for about 20 
seconds. The first 500 microseconds of this spih was "fast" extraction for neutrino 
experiments. A typical run lasted 40 minutes. 
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Fig. 14. Rate dependences of several quantities versus the scaler deduced intensity 
monitor within a single accelerator spill. Figure (a) shows the average number of 
embedded pair tracks versus the beam intensity. Figures (b), (c), and (d) show the 
accidental rate in three Cerenkov cells versus the beam intensity. 

FOCUS used three multicell Cerenkov counters, operating in threshold mode. 
We believe that this algorithm may prove useful in future experiments em- 
ploying threshold Cerenkov counters in high rate environments. This article 
describes the CITADL algorithm, illustrates its effectiveness on charm sig- 
nals, and discusses the method we used for continuous monitoring of each 
cell's photoelectron yield and noise. 

Although the pulse heights for all of the firing Cerenkov cells in a given event 
were recorded, our algorithm did not use this pulse height information. The 
CITADL algorithm was based on the probabihty that Cerenkov cells uniquely 
associated with a given track either fired or failed to fire. Because there were 
only two outcomes per Cerenkov cell, we say the algorithm returned a "digital 
likelihood" for a given particle hypothesis. Given the large number of cells 
with different accidental rates and photoelectron yields comprising the FOCUS 
Cerenkov system,the "digital" likelihood provided an essentially continuous 



27 



identification variable. 

By using only the on/off status of the Cerenkov cells, we found that it was 
possible to include the effects of "accidental" firing due (in our case) to the 
untracked electromagnetic debris in regions close to our photon beam. Es- 
sentially we combined the probability of a cell firing due to Cerenkov light, 
with an accidental firing probability using De Morgan's Law. The accidental 
probability for each Cerenkov cell was parameterized in terms of an intrinsic 
accidental rate and a contribution proportional to the instantaneous beam 
intensity. By realistically including the accidental rate in our likelihood, we 
substantially improved our ability to identify light particles (such as pions) 
over the less sophisticated algorithm used in our previous charm photoproduc- 
tion experiment, E687. Interestingly enough, the ability to positively identify 
pions with high efficiency and low kaon contamination, proved useful in signif- 
icantly increasing our signal to noise even in Cabibbo favored charm decays. 

The goal of the CITADL algorithm was to provide fiexible identification with a 
broad efficiency versus misidentification curve. The fiexibility of this algorithm 
has proven very useful in assessing systematic errors due to misidentified charm 
refiections in our recent studies of Cabibbo suppressed and doubly suppressed 
charm decays. [4] -[6] 
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